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The low-energy theory for multi-wall carbon nan- 
otubes including the long-ranged Coulomb interactions, in- 
ternal screening effects, and single-electron hopping between 
graphite shells is derived and analyzed by bosonization meth- 
ods. Characteristic Luttinger liquid power laws are found for 
the tunneling density of states, with exponents approaching 
their Fermi liquid value only very slowly as the number of 
conducting shells increases. With minor modifications, the 
same conclusions apply to transport in ropes of single-wall 
nanotubes. 

PACS numbers: 71-10.Pm, 71.20-Tx, 72.80. Rj 



Metallic carbon nanotubes constitute a novel and ex- 
citing realization of one-dimensional (ID) conductors 
plpl . The strong electronic correlations observed experi- 
mentally [pf pose many challenges to theorists. It is well 
known that electron-electron interactions invalidate the 
Fermi liquid description in one dimension. Often ID con- 
ductors can instead be described as a Luttinger liquid 
(LL) at low energy scales. The theoretical prediction [|||] 
of LL behavior in a single- wall nanotube (SWNT) has in- 
deed been verified in a recent transport experiment J3|. 
Furthermore, pronounced interaction effects have been 
observed in experiments on individual multi-wall nan- 
otubes (MWNT) composed of several (typically 10 to 
20) concentric graphite shells Ma, most notably a pro- 
nounced zero-bias anomaly at low applied voltage, which 
was tentatively interpreted in terms of LL theory. Since 
MWNTs are much easier to manipulate than SWNTs, 
interaction effects might be useful for applications, and 
whether a MWNT can display LL behavior or not could 
then be a question of practical importance. Obviously, 
very similar questions [and answers given below] apply to 
bundles (ropes) of SWNTs, the system actually studied 
in Ref. §. 

In this Letter, the low-energy theory of an individual 
MWNT composed of N metallic graphite shells with 
radii R\ < R2 < ■ ■ ■ < Rn is derived, fully taking into 
account the externally unscreened Coulomb interaction, 
internal screening effects, and inter-shell electron tunnel- 
ing. Insulating shells and the substrate are incorporated 
in terms of a (space-dependent) dielectric constant. The 
theory holds for energy scales fcsT, eV -C v/ Rn [we set 
% = 1] , where v sw 8 x 10 5 m/sec is the radius- independent 
Fermi velocity; for Rn — 10 nm, v/Rn ~ 350 meV. Then 
one needs to take into account only two transport bands 
a = ± per shell. For both bands, there is a right- and a 



left-moving (r = ± = R/L) branch with linear (massless) 
dispersion. Regarding the importance of disorder, both 
theory pi] and experiment |3MO are consistent with 
large elastic mean free paths (several hundred nm up to a 
few fiva) in metallic tubes. Therefore we mainly focus on 
the clean case here. By employing bosonization methods 
Jl2| |, pronounced LL effects are predicted for a MWNT, 
with only a slow crossover to Fermi liquid behavior as N 
increases. The precise relation between experimentally 
measurable exponents and microscopic quantities turns 
out to be quite different from the SWNT case discussed 
in Refs. |,§. 

The construction of the low-energy theory starts by 
expanding the electron operator for spin a — ± on shell 
n = 1, . . . , N in terms of the right- and left-moving Bloch 
waves (j)ran{x,y) for each transport band a = ±, 
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The Bloch waves are stiff around the circumference, 
4>ran(x,y) = (2nR n )~ 1 ' 2 exp(— iakpx) , where kp — 
2ir/3a with the lattice constant a = 2.46 A, < y < 
2irR n , and x is the transport direction. The expansion 
(0) allows to formulate the theory in terms of ID fermion 
operators i/Vao-n {x) . Ignoring interactions and inter-shell 
tunneling for the moment, the linear dispersion implies 
N copies of a ID massless Dirac Hamiltonian, 
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In the next step, let us include the Coulomb interac- 
tions among the electrons. Short-ranged backscattering 
processes are ignored since the corresponding couplings 
are extremely small even in a SWNT and scale as l/R n 
H . Since one is normally off half- filling due to the pres- 
ence of metallic gates, Umklapp scattering is also ne- 
glected. Under the expansion (|l]) , the important forward 
scattering processes lead to the contribution 

1 N f 
-ff / = 2 E / dxdx' p n (x)V nm (x- x')p m (x') , (3) 

n,m— 1 

where the ID densities for shell n are 

pn(x) = Y^ ^i aa ni x )^r a <jn{ x ) ■ ( 4 ) 

raa 

The effective ID potential V nm = V mn is obtained from 
the externally unscreened 3D Coulomb potential by in- 
tegrating over the circumferential coordinates using the 
Bloch functions, 
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The dielectric constants n nm include the effect of the sub- 
strate and of insulating shells. The Fourier transform 
V nm (k) at long wavelengths |fcitV| <C 1 takes the form 

2e 2 
V nm (k) = -—\n(\k\R nm ), (5) 

Knm 

where the "mean radius" R nm reads 

\nR nm = -r ln{[(i?„ - R m )/2} 2 + R n R m sin 2 V } . 

Jo 4tt 

An order-of-magnitude estimate can be obtained by re- 
placing sin 2 ip -► 1/2, yielding R„ rn « [(i? 2 + R 2 .)^] 1 / 2 . 
The logarithmic singularity of V nm (k — > 0) reflects the 
unscreened 1/r tail of the Coulomb interaction. In order 
to make contact with the usual LL concept, we restrict 
ourselves to a long-wavelength description by cutting off 
this singularity at k = 2tt/L, where L is the MWNT 
length |l5| ]. Including a factor 4/ttv for later convenience, 
this leads to the dimcnsionless coupling constants 



Un 



T^VKjuyi 



HL/2irR nm ) 



(6) 



For L = \ /ira, R nm = 6 nm and K„ m = 1.4, Eq. (J6J) yields 
U nm ~ 16. Therefore both the intra- and the inter-shell 
electrostatic interactions are typically very strong and 
of comparable magnitude. Due to the weak logarithmic 
dependence of U nm on R nm , the approximation U nm = U 
can already yield sensible results. 

By virtue of bosonization, Hq + Hi can now be diago- 
nalized. The bosonized form of the ID fermion operators 

ipra<rn(x) = '^- expi iq F rx + i(7r/4) 1/2 [0 c+! „ + rip c+ , n 
y2na ( - 

+ uQ c -. n + ra</J c _,„ + (j9 s+tn + ra(p s+t „ 

+ aa9 s -. n + raaip s -^ n } j , (7) 

where symmetric and antisymmetric linear combinations 
of the two transport bands a = ± have been formed for 
charge and spin degrees of freedom. The resulting four 
channels are labelled by the index 7 = (c+, c— , ,s+, s— ). 
The boson phase fields obey the commutator algebra 



[0^ in (x),ipy t n'(x')] = -(i/2)5 nn >5~ /r 
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so that ip-y y n has the canonical momentum il 7i „ = 
—d x ly „. The Majorana fermions r] raan ensure anticom- 
mutation relations between ID fermions with different in- 
dices raan. Finally, qp — Ep/v is determined by exter- 
nal gate voltages. Due to inter-shell tunneling, the Fermi 



energy Ep and hence qf are identical for all N shells. 
Using Eq. (0), the density is p n (x) = (4/7r) 1 / 2 9 x ^ c+i „, 
and as a consequence the Hamiltonian decouples in all 
four channels, Hq + Hi = J2y H 1 . The charged channel 
is described by 
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The three neutral channels correspond to Eq. (||) with 
U nm — 0. Diagonalizing H c+ then leads to an eigen- 
value problem similar to the one studied by Matveev and 
Glazman [14] for many-channel quantum wires in semi- 
conductor heterostructures, 
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where (p c + >n (x) = ^ ■ T n j$j(x) and T n j is an orthogonal 
matrix. In terms of the new boson phase fields $.,■ and 
their canonically conjugate momenta lij, H c+ results in 
the standard LL form, 



H r 



N r 
J2 dx{u 2 +gr 2 (d x ^y 



(10) 



The LL interaction constants g,j < 1 for the N eigen- 
modes described by the phase fields $_,■ serve as a measure 
of the correlation strength [021 . 

It is then easy to determine all scaling exponents of 
interest. The exponents ri b / e of the tunneling density of 
states (TDOS), p(E) ~ E' n , for tunneling of an electron 
into the outermost shell [n — N) in the bulk or close to 
the end of the tube are 
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These exponents govern the power laws G ~ T 11 of the 
temperatur-dependent linear conductance for tunneling 
into the MWNT and can be measured in the experimen- 
tal setup of Ref. J7|], where external leads contact only 
the outermost shell of the MWNT. Moreover, in the limit 
of weak disorder backscattering, the linear conductance 
corrections 5G ~ T~ p are characterized by the exponent 
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At sufficiently low temperatures, disorder in a LL always 
leads to a strong backscattering situation. Then the con- 
ductance vanishes as T° for T — v 0, where j3 = 2rj e . Ob- 
servation of the exponents p and (3 requires good contacts 
between external leads and the MWNT. 



It follows from the above analysis that an individual 
MWNT exhibits LL behavior on low energy scales, where 
the LL exponents can be computed by solving the eigen- 
value problem (|9|). Of course, exactly the same reason- 
ing applies to ropes of SWNTs, the only difference being 
the precise form of the couplings U nm . Indeed the log- 
arithmic singularity in Eq. (o) and hence the order of 
magnitude of the couplings applies to both systems. For 
N — 1, all expressions above reduce to the SWNT theory 
of Refs. |,§. 

To make analytical progress, we now solve Eq. (m for 
the special case where all U nm are equal and given by 
U. [In the estimates below, we neglect the weak loga- 
rithmic iV-dependence of Eq. (y) and put U — 16.] The 
eigenvalues are 

gi=g=(l + NU)- 1 ' 2 , 92,...,n = 1. (14) 

The first value has eigenvector r„i = iV -1 / 2 and is iden- 
tified with the standard collective LL plasmon mode, now 
characterized by a pronounced iV-dependence of the LL 
interaction parameter g. The remaining N — l degenerate 
modes are not affected by the interactions and correspond 
to Fermi liquid quasiparticles. From the orthogonality of 
T n j, it follows that J2j>i ^nj = 1 — V-^i an d the above 
exponents read 
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For N — > oo, the exponents rjb^ e and p approach zero (the 
Fermi liquid value) only as iV -1 / 2 , because the LL pa- 
rameter g also goes to zero as N~ x ' 2 . This slow approach 
implies that LL power laws should be well observable in 
typical MWNTs by focusing on the TDOS. For instance, 
for N — 10 the exponents are n e — 0.292 and r/i, = 0.135, 
while the corresponding SWNT values are n e — 0.780 and 
rjb = 0.295. On the other hand, the exponent p scales as 
1/JV for large N and will thus vanish more rapidly. For 
instance, for N = 10 it is p = 0.046, but for N = 1 we 
would get p — 0.379. LL power laws in MWNTs or ropes 
of SWNTs are thus generally much more pronounced for 
the TDOS than for the backscattering corrections in the 
presence of weak disorder. This important fact may be 
used in practice to obtain information about the number 
N of conducting shells. 

So far single-electron hopping between the shells has 
been ignored. It can be written as 



N 
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where the hopping matrix acts only between nearest- 
neighbor shells, T nm = -t n 6 n +i,m - t n -iS n ,m+i- The 



hybridization matrix elements are of the order t n « 
250 meV 0. Processes where raa-type fermions are 
scattered into different racr-states are suppressed against 
Eq. (nS) by momentum conservation and by higher scal- 
ing dimensions. Renormalization group arguments sug- 
gest that H t is relevant, so one has to be careful about 
its effect. 

To include the hopping, we go back and first diagonal- 
ize Hq + H t . With the orthogonal matrix Q nv , 
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the rotated fermions ^ ra au again obey the ID Dirac 
Hamiltonian (g). Denoting the N eigenvalues of T nm 
as T v , the respective eigenvectors span Q nv so that 
H t = ~Y^ V T V f dxp„(x), with ID densities p„ defined in 
analogy to Eq. (g). In the long- wavelength limit, Eq. (||) 
now reads 
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with matrix elements 



-^1/11/2^3 ^4 — / j^ nm^t/inWn-^^/^m^f 
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The complicated four-fermion interactions (EOh are to a 
large degree responsible for the technical difficulties en- 
countered in previous studies of coupled LLs ]12[] , Ac- 
cording to Eq. (ft), however, in MWNTs the couplings 
U nm are approximately equal, which allows to make fur- 
ther progress. The situation might be different in ropes 
of SWNTs, where one should first diagonalize Ho + Hj 
by solving Eq. (ph, and only afterwards include H t , e.g., 
by perturbation theory. 

In MWNTs, the dominant matrix elements in Eq. ( pl| ) 
are Y vvv > v > = W vv > = W v 'v In fact, if U nm = U, or- 
thogonality of Q nl , implies that W vv i = U and all other 
matrix elements in Eq. (|2l]) are zero. Keeping only the 
couplings W uv ' j Eq. d2G) simplifies to 
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(22) 



Remarkably, the model Ho + Hi + H t can then be solved 
exactly by bosonization in the rotated basis. Using 
Eq. (Q) for ipraav(x), the neutral channels are described 
by the same Hamiltonian as before. The only change 
arises in the charged sector, 



N 



dx (IT 2 + <?- 2 (d x %) 2 + V e j N j , (23) 
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where Cj — J2 u T u T„j. The "zero modes" Nj — 
(4/7T) 1 / 2 J dxd x <&j denote the total number of particles 
in the jth eigenmode. The eigenvalues gj and the ma- 
trix T v j are then determined from Eq. (0) with U nm be- 
ing replaced by W vv > . For U nm = U, the only yet im- 
portant effect of the hopping consists of a splitting of 
the iV ID transport bands, with all plasmon excitations 
and LL interaction constants gj being unaffected. Using 
J2j F^jej — T u , the TDOS for tunneling into the outer- 
most shell is generally of the form 
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where p v (E) ~ E nu > h / e . The bulk/end exponents 77„.b/ e 
for tunneling into the z/th MWNT eigenmode are given 
by Eqs. ( |il"| ) and (|l^), respectively, with the replacement 
^Nj —> r„j. For U nm = U, both exponents are indepen- 
dent of v and given by Eqs. ( |l5| ) and (|l6| ). 

Next the other interactions in Eq. (EQ) are briefly an- 
alyzed. If three or four indices Vi are different, the re- 
sulting perturbations are always highly irrelevant since 
single-particle excitations for more than two modes are 
involved in the scattering. The contribution with v\ = v-$ 
and Vi — 1/4 corresponds to a simultaneous transfer of 
two fermions from one mode to another, which is also an 
irrelevant process. The remaining interactions are due 



liquid thus proceeds by (i) a iV -1 / 2 decrease in the ex- 
ponents r]b/ e , and (ii) tunneling thresholds that become 
closer and closer spaced. For moderate values of N, how- 
ever, different T v are sufficiently well separated to allow 
for the observation of LL power laws. 

To conclude, the low-energy theory of an individ- 
ual MWNT on an insulating substrate has been given. 
Including the long-ranged Coulomb interaction, pro- 
nounced correlation effects are predicted which can be 
understood in the framework of the Luttinger liquid 
model. Clearly, there are many open questions to be 
addressed in future work, e.g., the effects of a parallel or 
perpendicular magnetic field, or the virtual or real pop- 
ulation of higher subbands as the energy scale increases. 

I wish to thank L. Forro, A. Gogolin, H. Grabert, and 
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support by the Deutsche Forschungsgemeinschaft (Bonn) 
under the Gerhard-Hess program. 



to the couplings B vv i = B v 
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which are small 



for Unm ~ U. Using bosonization, one observes that the 
B vv i are marginally relevant scaling fields causing a weak 
renormalization of the LL interaction parameters also in 
the neutral channels [which were all equal to one before] , 
and the nonlinear contribution 



H'x = 
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Since marginally relevant scaling fields lead to exponen- 
tially small gaps for small couplings, the generic behavior 
of MWNTs is indeed described in terms of the LL model 
outlined above. The predicted pseudogap in ropes ]lq ] 
can then be rationalized by noting that the B vv i are pre- 
sumably larger in ropes. 

How will Fermi liquid behavior be restored as the 
number of conducting shells approaches infinity? Tak- 
ing for simplicity all nearest-neighbor hybridization ma- 
trix elements t n — t, the eigenvalues of the hopping 
matrix for large N read T„ = —2tcos(2TTi//N), with 
QnvQn„ = l/N. Hence the TDOS becomes 



p(E)^N- 1 J2(E-T 1/ ] 
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where the threshold energies T v become narrowly spaced. 
In effect, for large enough N, the LL singularities will 
then be smeared out simply because subsequent thresh- 
olds are too close together. The approach to the Fermi 
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